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New Measurements of Orbital Period Change in Cygnus X-3 
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Abstract. A nonlinear nature of the binary ephemeris of Cygnus X-3 indicates either a change in the orbital 
period or an apsidal motion of the orbit. We have made extended observations of Cygnus X-3 with the Pointed 
Proportional Counters (PPCs) of the Indian X-ray Astronomy Experiment (IXAE) during 1999 July 3—13 and 
October 11—14. Using the data from these observations and the archival data from ROSAT, ASCA, BeppoSAX 
and RXTE, we have extended the data base for this source. Adding these new arrival time measurements to 
the published results, we make a comparison between the various possibilities, (a) orbital decay due to mass loss 
from the system, (b) mass transfer between the stars, and (c) apsidal motion of the orbit due to gravitational 
interaction between the two components. Orbital decay due to mass loss from the companion star seems to be 
the most probable scenario. 
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1. Introduction 

The nature of the compact object in the bright X-ray bi- 
nary Cygnus X-3, is a subject of much debate. In spite of 
being one of the most frequently observed X-ray sources, 
presence of an X-ray pulsar, a black hole, or a low mag- 
netic field neutron star has not yet been established. There 
is also uncertainty about the mass and type of the com- 
panion star. An interesting way to probe this system is 
to investigate the arrival time history of the 4.8 hr orbital 
modulation in the X-ray light curve. 

The unusual X-ray binary Cygnus X-3 is located in the 
plane of our galaxy at a distance of > 11.6 kpc (Dickey 
1983). Because of the high luminosity of the source in X- 
ray, infrared and radio bands, it has been observed on 
many occasions at these wavelengths. Strong optical ex- 
tinction in the direction of the source prevents optical ob- 
servations. The mass of the companion star has been de- 
termined by many different means and it is found to be in 
the range of a fraction of solar mass to a few solar mass 
(Van den Heuvel & De Loore 1973; Tavani, Ruderman, & 
Shaham 1989; van Kerkwijk et al. 1992). Parsignault et al. 
(1976) discovered a periodicity of 4.8 hr in the X-ray flux 
with nearly sinusoidal variation, which is believed to be 
due to the orbital motion. The variation in the sinusoidal 
shape of the light curve from cycle to cycle is reported by 
van der Klis and Bonnet-Bidaud (1981). They showed that 
the shape of the average light curve formed from succes- 
sive cycles is quasi-sinusoidal with a slow rise and fast fall. 
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Several models have been proposed to explain the deep 
and near-sinusoidal modulation of the X-ray intensity of 
Cygnus X-3 with the orbital period (Ghosh et al. 1981). 
Material in the form of a large shell around the binary sys- 
tem can produce the X-ray modulation by scattering, ab- 
sorbing, and re-emitting the X-rays from the compact ob- 
ject, if the X-ray emission is only from the non-shadowed 
part of the shell. White & Holt (1982) proposed the ac- 
cretion disk corona (ADC) model to explain the modu- 
lation in the X-ray light curve of the source. According 
to this model, the X-ray source is covered by an optically 
thick corona with a radius of about 10 9 cm, an optical 
depth of about 10 and a temperature of about 2 keV. 
Another model which tries to explain the orbital modu- 
lation is stellar wind model (Willingale, King, & Pounds 
1985; Kitamoto et al. 1987). This model assumes a strong 
and highly ionized stellar wind from the companion star 
with an optical depth of about 1 and a radius of about 10 11 
cm and the temperature of the wind is model dependent. 

Evolution of the orbit of Cygnus X-3 is studied by 
measuring the arrival times of the minima in each orbital 
motion. The ephemeris of 4.8 hr modulation in the X-ray 
light curve of the source has been studied by many authors 
(Leach et al. 1975; Mason & Sanford 1979; Parsignault et 
al. 1976; Lamb, Dower & Fickle 1979; Eisner et al. 1980; 
van der Klis & Bonnet-Bidaud 1981, 1989; Kitamoto et 
al. 1987). The time derivative of the orbital modulation 
period (P) and second derivative (P) have been measured 
as ~ 10~ 9 s s _1 and ~ -10 -11 yr (van der Klis & Bonnet- 
Bidaud 1981, 1989; Kitamoto et al. 1987). The unusually 
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large value of P/P (= 2.2 x 10~ 6 yr" 1 ), which if linked 
to the binary evolution of the system, is similar to several 
short-period bright LMXBs such as X1822-371 with P/P 
~ 3.4 x 10~ 7 yr _1 and orbital period of 5.57 hr (Hellicr 
et al. 1990; Parmar et al. 2000), EX00748-767 with P/P 
~ 2 x 10~ 7 yr _1 and orbital period of 3.82 hr (Parmar et 
al. 1991). This value of P can be interpreted as due to the 
mass transfer from the optical companion to the compact 
object. Using the measured period derivative of Cygnus 
X— 3, Kitamoto et al. (1987) determined the rate of mass 
loss from the binary companion as about 10~ 6 M Q yr _1 . 
Tavani et al. (1989) tried to explain the observed high 
value of P/P in Cygnus X-3 and the required mass trans- 
fer rate close to — 1.58 x 10~ 8 M© yr _1 . They claimed 
that the companion is a degenerate star with solar com- 
position and in the mass range of 0.01 < m < 0.03 M Q 
which under-fills its Roche lobe. The X-ray illumination 
from the primary compact object produces the evapora- 
tive wind from such a low-mass degenerate companion. 
However, from infrared observations, van Kerkwijk et al. 
(1992) discovered that the binary companion of Cygnus X- 
3 has the spectrum of a Wolf-Rayet star and predicted the 
mass of the companion as ~ 10 M Q which suggests that 
Cygnus X-3 is a high mass X-ray binary. Recent obser- 
vations by Fender et al. (1999) confirm the nature of the 
companion as an early-type WN Wolf-Rayet star. 

We have measured the available arrival time of min- 
imum value of the light curve from two observation of 
Cygnus X-3 with PPCs of IXAE in 1999 and several 
archival data. By combing these new measurements with 
the previously published results, we are able to determine 
the rate of change of orbital period with greater accuracy. 

2. Observations and Archival Data 

2.1. PPC Observations 

The X-ray observations of Cygnus X-3 were carried out 
twice on 1999 July 3-13 and 1999 October 11-14 with 
1 s time-integration mode with the PPCs of the IXAE 
on board the Indian satellite IRS-P3. The IXAE includes 
three co-aligned and identical, multi-wire, multi-layer pro- 
portional counters with a total effective area of 1200 cm 2 
covering 2 to 18 keV energy range with an average de- 
tection efficiency of about 60% in 3—7 keV energy range. 
Each PPC is equipped with a honeycomb type collimator 
having a field of view of 2°. 3 x 2°. 3 and has an overall 
energy resolution of about 22% at 6 keV (Agrawal et al. 
1998; Rao et al. 1998). The IRS-P3 satellite is in a circu- 
lar orbit at an altitude of 830 km and inclination of 98° . 
Pointing towards any particular source is done by using a 
star tracker with an accuracy of < 0°.l. The useful obser- 
vation time is limited to the latitude range typically from 
30°S to 50°N to avoid the rise in charged particle back- 
ground at high latitudes. The South Atlantic Anomaly 
(SAA) region restricts the observation to 5 of the 14 or- 
bits per day. As the detectors are co-aligned, simultaneous 
background observation is not possible. The background 



counts are measured after the observation of a source by 
pointing the PPCs to a source-free region in the sky, close 
to the target source. The total useful exposure time of 
the X-ray source during the two observations is 31245 s. 
The PPC light curves did not show strong intensity varia- 
tions on the time scale of a few seconds to a few hundred 
seconds. 

2.2. Archival data 

We have also analyzed archival data from several satel- 
lites. From ROSAT mission, we have analyzed HRI data 
with two different observation IDs. The durations of these 
two observations are inter-wined, but the overall count 
rate is different because in one of the observations, the 
source was near the edge of the field of view where effec- 
tive area is small. RXTE/PCA observations of the source 
in 2—60 keV energy range, combined data from GIS detec- 
tors (GIS2 and GIS3) and SIS detectors (SISO and SIS1) 
of the ASCA satellite and data from MECS detectors of 
BeppoSAX observations were used to determine more ar- 
rival time information of Cygnus X— 3. We have analyzed 
the RXTE/ASM dwell data in the energy range of 2-12 
keV for the entire period of observation starting from 1996 
January to 2001 August 30. The details of the observations 
of the source with different satellites are given in Table 1 
along with the useful exposure period. 

Table 2. Archival data used for the determination of the 
arrival time of the minima of the binary modulation 



Satellite Date of observation Exposure Detector 

(in seconds) 



ASCA 


1994 


29/5 - 


30/5 


21468 


GIS 


ASCA 


1994 


29/5 - 


30/5 


17527 


SIS 


ROSAT 


1995 


20/4 - 


22/4 


72709 


HRI 


ROSAT 


1995 


21/4 - 


23/4 


51670 


HRI 


BeppoSAX 


1996 


22/9 - 


24/9 


51626 


MECS 


BeppoSAX 


1996 


10/10 




26167 


MECS 


BeppoSAX 


1996 


29/10 




13664 


MECS 


BeppoSAX 


1997 


19/9 - 


20/9 


25802 


MECS 


RXTE 


1996 


24/8 - 


30/8 


86450 


PCA 


RXTE 


1997 


16/2 - 


22/2 


59465 


PCA 


RXTE 


1997 


05/6 - 


26/9 


64693 


PCA 


RXTE 


1998 


16/5 - 


21/5 


34723 


PCA 


RXTE 


1996 
2001 


04/1 - 
30/8 






ASM 



3. Data Analysis and Results 

3.1. Arrival time determination 

The data used for the analysis were corrected to the solar 
system barycenter. To avoid the cycle to cycle variation in 
the light curve, we have created the time averaged orbital 
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Fig. 1. Folded light curves in 2—18 keV and 6—18 keV energy ranges are shown in the panels (a) and (b) (for 1999, 
July 3—13 PPC observations) and panels (d) and (c) (for 1999, October 11 — 14 PPC observations) on a normalized 
scale with the template data (described in text) which is indicated by lines. The panels (c) and (f) show the folded 
hardness ratios (ratio of the count rate in 6—18 keV energy range to the count rate in 2—6 keV energy range for 1999, 
July and October observations respectively. 



modulation profile by folding the source light curves with 
a period of 17253 seconds (Eisner et al. 1980). In all cases, 
the shape of the folded light curves is somewhat similar 
to the shape of the template function which is the long 
term average shape of the light curve obtained with the 
Copernicus satellite (Mason et al. 1976). 

To determine the arrival times of the minimum of the 
folded light curves, we have used the normalized template 
data of X-ray ephemeris of Cygnus X-3 (van der Klis et 
al. 1989). The folded data were cross-correlated with the 
standard template function as described by van der Klis 
et al. (1989) to get the phase difference in the minima of 
the template and the data which give the measure of the 
arrival time of the minima. In addition to the statistical 
variations, Cygnus X-3 often shows some short-term in- 
tensity variations superposed on the orbital modulation. 
We have multiplied the statistical error estimates of each 
folded light curve by some constant factors to make them 
fit to the template with a reduced x 2 of 1.0. The error 
in the arrival times were estimated from the uncertain- 
ties in determining the center of the Gaussian function of 
the cross-correlated data. The short-term occasional inten- 
sity variations also cause some random fluctuations in the 
arrival time determinations. Following the earlier works 
(Kitamoto et al. 1995), we have quadratically added 0.002 



d to the arrival uncertainties to make them compatible 
with the previous reported measurements. 

This analysis yielded two new arrival points from PPC 
observations (Singh et al. 2001), two points from ASCA, 
two points from ROSAT, five points from RXTE/PCA, 
and four points from BeppoSAX observations. We divided 
the RXTE/ASM light curve into 19 segments of durations 
of 100 days as given in Table 2, and determined the ar- 
rival times from each data segment. The information re- 
garding the new arrival times are given in Table 1 along 
with the observation periods of the source with different 
instruments on different satellites and the orbit number. 
In accordance with the earlier works (Leach et al. 1975), 
we have assigned the orbit number of the source as at 
JD 2440949.9201. The orbit numbers for later observa- 
tions were calculated by using 4.8 hour orbital period of 
the source. The significant increase in the number of the 
new arrival points to the old data improved the numerical 
values of the period derivative and the double derivative. 

The orbital modulation of the light curves obtained 
from the observations of the source with different instru- 
ments are shown in different figures along with the nor- 
malized template. Panel (a) and (d) of Figure 1 show the 
folded light curves in 2—18 keV energy range along with 
the normalized template for 1999 July 3—13 and 1999 
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Phase 

Fig. 2. The figure shows the normalized folded light curves 
of the source obtained from the RXTE/ASM data for dif- 
ferent energy ranges and at different time. Upper pan- 
els (panels a, b, and c) show the folded light curves in 
the energy ranges 1.3—3.0 kcV, 3.0—5.0 keV, and 5.0-12.1 
kcV for the data in the time range JD 2450086.8222 - 
2450186.4878 whereas the lower panels (panels d, e, and f) 
show the folded light curves in the energy ranges 1.3—3.0 
keV, 3.0-5.0 keV, and 5.0-12.1 keV for the data in the 
time range JD 2450887.2862 - JD 2450986.1345 respec- 
tively. 

October 11 — 14 PPC observations of the source respec- 
tively. In Figure 2, we have shown the orbital modulation 
of the light curves along with the normalized template for 
two data segments of the RXTE/ASM for different energy 
bands. From these figures, it is observed that the shape 
of the orbital modulated light curves of Cygnus X-3, ob- 
tained from various observations, is similar to the shape 
of the normalized template. 

3.2. Orbital evolution of the binary system 

Figure 3 (top panel) reports the new measurements of the 
arrival time of the minima in the orbital modulation of the 
X-ray light curves of Cygnus X-3 together with those of all 
the previously published values (Kitamoto et al. 1987; van 
der Klis & Bonnet-Bidaud 1989). We have tried to fit the 
arrival time data with different polynomial functions. The 
residuals of the various polynomial fits performed on these 
data are also shown in other panels of the figure. The sec- 
ond panel of Figure 3 shows the residual of the linear fit to 
the arrival data whereas the third, fourth, and fifth panels 
of Figure 3 show the residuals of the quadratic, cubic, and 
sinusoidal fit to the arrival time data respectively. 

The results of the various polynomial fits performed on 
the complete set of arrival time data are listed in Table 
3. The arrival times of the orbital modulations clearly 
deviate from a linear relationship with the orbit num- 
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Fig. 3. The figure shows the arrival time history for the 
source Cygnus X-3 (top panel). The data points beyond 
orbit number 42830 are obtained from the present work. 
Arrival time residuals of the orbital modulation of the 
source with respect to the best-fit linear, quadratic, cubic 
and a linear and sinusoidal models to the data are shown 
in different panels. 



bers. Addition of a quadratic term fits the data much 
better with a reduced \ 2 of 3.08 (113 dof). From the 
quadratic fit, the value of the period derivative is derived 



-10 



i.e, an evolution time scale 
1 . These values are similar to 



to be P = 5.76333 x 10 
of P/P=l. 05419 x 10~ 6 yr" 
what was obtained with the data prior to 1993 (Kitamoto 
et al. 1995). A cubic model for the arrival time fits the 
data slightly better with reduced x 2 of 2.96 (112 dof). 
However, the rate of change of P is found to be much 
smaller compared to that obtained by Kitamoto et al. 



(1995), P =-1.34008 x 10" 



The curved feature of 



the residual of the linear fit to the data prompted us to fit 
with a model consisting of linear and sinusoidal terms as 
model components. Although the reduced \ 2 f° r this fit- 
ting is comparable to the reduced % 2 for cubic model, the 
amplitude of sinusoidal component is unphysical because, 
according to the method described by Batten (1973), it 
requires the eccentricity of the binary system to be >• 1. 
In the linear plus sinusoidal model, a reduced \ 2 of 3.5 
(dof 112) was obtained by restricting the amplitude of 
the sinusoidal component within a physical limit so that 
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Observation 


Satellite 


Number of 


Arrival time 


Error in 


Period 




(detector) 


Orbits 


(2440000 JD+) 


arrival time 


(MJD) 






(p= 17253 s) 


(day) 


(day) 


49501.91- 


49502.84 


ASCA(GIS) 


42832 


9502.8736 


0.0021 


49501.91- 


-49502.84 


ASCA(SIS) 


42832 


9502.8736 


0.0028 


49827.99- 


-49829.44 


ROSAT(PSPC) 


44466 


9829.1748 


0.0020 


49828.59- 


-49830.56 


ROSAT(PSPC) 


44470 


9829.9766 


0.0021 


50086.32- 


-50185.98 


RXTE(ASM) 


46006 


10136.6934 


0.0020 


50186.06- 


-50285.99 


RXTE(ASM) 


46506 


10236.5381 


0.0020 


50319.46- 


-50325.83 


RXTE(PCA) 


46940 


10323.2031 


0.0024 


50348.92- 


-50350.18 


BeppoSAX(MECSl) 


47074 


10349.9609 


0.0020 


50385.00- 


-50385.52 


BeppoSAX(MECS) 


47253 


10385.7090 


0.0025 


50386.31- 


-50485.96 


RXTE(ASM) 


47508 


10436.6299 


0.0020 


50495.02- 


-50501.12 


RXTE(PCA) 


47818 


10498.5391 


0.0047 


50486.02- 


-50585.87 


RXTE(ASM) 


48008 


10536.4727 


0.0020 


50586.16- 


-50685.95 


RXTE(ASM) 


48509 


10636.5234 


0.0020 


50604.75- 


50717.50 


RXTE(PCA) 


48635 


10661.6787 


0.0021 


50710.65- 


50711.26 


BeppoSAX(MECS2) 


48884 


10711.4014 


0.0020 


50731.36- 


-50731.90 


BeppoSAX(MECS2) 


48988 


10732.1738 


0.0021 


50686.01- 


-50785.73 


RXTE(ASM) 


49009 


10736.3643 


0.0020 


50786.51- 


-50883.86 


RXTE(ASM) 


49506 


10835.6104 


0.0021 


50886.79- 


-50985.63 


RXTE(ASM) 


50012 


10936.6533 


0.0020 


50949.63- 


-50954.98 


RXTE(PCA) 


50093 


10952.8203 


0.0035 


50986.41- 


-51085.89 


RXTE(ASM) 


50513 


11036.6934 


0.0020 


51086.09- 


-51185.95 


RXTE(ASM) 


51013 


11136.5381 


0.0047 


51186.22- 


-51285.96 


RXTE(ASM) 


51514 


11236.5830 


0.0029 


51286.03- 


-51385.62 


RXTE(ASM) 


52013 


11336.2324 


0.0021 


51361.49- 


-51371.77 


IXAE(PPC) 


52168 


11367.1836 


0.0021 


51404.44- 


-51404.64 


RXTE(PCA) 


52358 


11405.1309 


0.0022 


51386.16- 


-51485.96 


RXTE(ASM) 


52515 


11436.4844 


0.0023 


51461.72- 


-51464.79 


IXAE(PPC) 


52651 


11463.6338 


0.0024 


51549.73- 


-51649.73 


RXTE(ASM) 


53334 


11600.028320 


0.0020 


51649.90- 


-51749.80 


RXTE(ASM) 


53836 


11700.274410 


0.0020 


51750.06- 


-51850.02 


RXTE(ASM) 


54338 


11800.516600 


0.0020 


51850.08- 


-51948.81 


RXTE(ASM) 


54836 


11899.958010 


0.0020 


51951.46- 


-52051.45 


RXTE(ASM) 


55347 


12002.000980 


0.0020 


52051.58- 


-52151.53 


RXTE(ASM) 


55848 


12102.048830 


0.0020 



e < 1. However, if a precessing eccentric orbit is the rea- 
son for the non- linear nature of the arrival times, a change 
in the shape of the orbital modulation of the light curve 
is expected, which has not been observed in the last ~ 30 
years. Based on this argument, apsidal motion seems not 
to be an appropriate description of the arrival time data. 

4. Spectral variations with orbital phases in 
Cygnus X-3 

To investigate the spectral variations of the source during 
different orbital phases, we have analyzed the hardness ra- 
tios from the PPC observations with IXAE on 1999, July 
3-13 and 1999, October 11 - 14 and the RXTE/ASM 
data since the beginning of the observation. From the cal- 
ibration of the PPCs using Crab, it is found that the spec- 
tral data from PPC— 3 is more reliable. Hence to obtain 
the hardness ratio (ratio of the count rate in 6 — 18 keV 
energy range to the count rate in 2 — 6 keV energy range) , 
we have used the data from PPC— 3. To study the spec- 



tral behavior of the source at different orbital phases, we 
have created the time averaged profile by folding the light 
curves and the hardness ratio data with a period of 17253 
seconds. We have plotted, in Figure 1, the orbital mod- 
ulated light curves in 2 — 18 keV and 6 — 18 keV en- 
ergy ranges and the hardness ratios for the 1999, July and 
October observations of the source. 

The folded RXTE/ASM light curves in 3-5 keV and 
5—12.1 keV energy bands plotted on a normalized scale 
with the template data and the corresponding hardness 
ratio for two different epochs (low and high states) are 
shown in Figure 3. From the PPC observations (Figure 
1) and the RXTE/ASM observations (Figure 4), it is ob- 
served that, a strong binary phase dependence of the spec- 
trum is noticed. 

5. Discussion 

The nature of the compact object in the binary system of 
Cygnus X-3 is not yet clearly understood. There are var- 
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Table 4. Statistics of different fitting models to the data of arrival time 



T ■ l ■ a r\r\r\ c -i -i A l C 

Linear ephemeris: \ = 4099 tor 114 dot 




T n = To + Pn 




To = 2440949.863 ± 0.002 JD 




T) a ififi/?077t; i a finnnnnn'7 J 

P = 0.19968775 ± 0.00000007 d 


Cov(io, r) = -1.6 X 10 


Parabolic ephemeris: \ 2 = 348.3 for 113 dof 




T n = To + P n + cn 2 


Cov(To, Po) = —2.3 x 10 a 


where c =P P/2 


Cov(T ,c) = 3.5 x 10 _16 d 2 


T) = 2440949.892 ± 0.001 JD 


Cov(Po, c) = —3.8 x 10 d 


Po = 0.19968443 ± 0.00000009 d 


P = (5.76 ± 0.24) xl0~ 1() 


c = (5.75 ± 0.16) x 10 d 


T) 1 7~> / 1 ACT _l_ A A/lA v/1A — 6„ „ — 1 

F 1 Po — (1.05 ± 0.04) xlO yr 


Cubic ephemeris: \ 2 = 331 for 112 dof 




T n = T + Ton + con + an 


s~\ / rri J~j \ ( Al\«,1A — 11 j2 

Uov(Jo, "o) = (—9.1) x 10 a 


where Co =PoP '/2 and 


Cov(To,c ) = 3.1 x 10~ 15 d 2 


d ~ To P/6 


Cov(T ,d) = (-3.2) x 10~ 20 d 2 


T = 2440949.8944 ± 0.0015 JD 


Cov(P ,c ) = (-5.5) x 10~ 19 d 2 


P = 0.199684009 ± 0.0000003 d 


Cov(P ,d) = 5.9 x 10~ 24 d 2 


rr, — (774-1 1\ v 1 n _11 W 

CO — 1 l.( I l.lt A 1U Ub 


Cnvlr,, rl\ — < 9 A\ X ~\(\~ 2S rl 2 
V_^OV^C0,tiJ — ^ ^•^) x 1U Ui 


d= (-2.4 ±1.3) x 10~ 16 d 


P = (7.7 ± 1.5) x 10- 10 




P= (-1.3 ±1.6) x lO -11 ^- 1 


Linear + Sinusoidal ephemeris: \ 2 ~ 391.1 for 112 dof 




T n = To + P n + c sin(2ir(n - 4>)/P aps ) 




T =2440949.92 ± 0.03 JD 




To = 0.199687 ± 9 x 10^ 07 d 




O = 0.059 ± 0.013d 




<f> = 334500 ± 6000 no. of orbit 




P aps =75.4 ± 0.2 yr 





ious arguments in the literature supporting the compact 
object as a massive black hole whereas some authors be- 
lieve the primary companion as a neutron star. Schmutz 
et al. (1996) suggested that the observed time variations 
in the profile of infrared emission lines from the binary 
system are due to the orbital motion of the companion 
which is a Wolf-Rayet star and derived mass function for 
Cygnus X-3 as 2.3 M Q . They obtained a range of 7-40 M 
as the mass of the compact object with a most likely value 
of 17 M Q and suggested that the binary system Cygnus X- 
3 is composed of a black hole (BH) and a Wolf-Rayet (WR) 
system. Several models have also been proposed to explain 
the lengthening of the orbital period (Bonnet-Bidaud & 
Chardin, 1988). Conservative mass transfer from the com- 
panion star can not explain the observations. This is due 
to the fact that mass accretion rate needed to produce 
the amount of energy emitted in X-rays and the observed 
orbital period change rate requires a very low mass com- 
panion ^0.02 M©, which will under- fill its Roche lobe to 
continue mass transfer. A low mass companion star with 
a moderate non-conservative mass loss rate, however, can 
explain the period change rate, if the mass loss is due to 
the X-ray heating of the atmosphere. On the other hand, if 
the companion star is massive, the orbit may widen only 
if there is significant mass loss from the system in the 
form of strong wind and the out-flowing mass carries spe- 
cific angular momentum of the primary star. The emission 
line rich X-ray spectrum of Cygnus X— 3 detected with 
Chandra (Paerels et al. 2000) indicates the presence of 



strongly photo-ionized stellar wind reminiscent of a high 
mass companion. 

6. Conclusion 

— We have measured 34 new minima and extended the 
measurement base to 30 years. 

— The parabolic model provides the best fit to the new 
data with a modified value of P. 

— Sinusoidal model and the corresponding apsidal mo- 
tion scenario is ruled out. 

— Mass loss from a massive companion is most favoured 
in Cygnus X— 3. 
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